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ON THE INFERENCE OF THE PHYSICAL PROPERTIES 
OF THE JOVIAN ATMOSPHERE FROM PHOTOMETRY 
OF ECLIPSES OF THE GALILEAN SATELLITES* 


by 

Michael J. Price 
Goddard Space Flight Center 


CHAPTER I 
INTRODUCTION 

Over 20 years ago Kuiper (1947) , reporting visual observations of the integrated 
brightnesses of Io and Europa as they entered eclipse, called attention to the existence of 
a tail on each light curve near the 14th magnitude. Kuiper tentatively attributed these 
tails to refraction in the Jovian atmosphere. This conclusion prompted the first ac- 
curate photoelectric photometry of the eclipses of the Galilean satellites, which were 
made by Kuiper in 1950 and by Kuiper and Harris in 1954. These observations and their 
interpretation have been discussed in detail by Harris (1961). 

From their observations , Kuiper and Harris attempted to determine the refraction 
properties of the Jovian atmosphere to infer its physical structure. Severe difficulties 
were encountered, however, in the interpretation of the data. Uncertainties in the radii 
of the satellites and in the distribution of reflectivity over their disks made the analysis 
arbitrary. Largely because of these difficulties, no further attempts have been made to 
use eclipses of the Galilean satellites to probe the Jovian atmosphere. 

In this paper the problem of inferring the physical properties of the Jovian atmos- 
phere from photometric observations of the eclipses of the Galilean satellites is thor- 
oughly discussed. All the theoretical and observational problems involved are examined 
in detail. It is shown that the difficulties encountered by Kuiper and Harris disappear if 
the photometry is extended to lower light levels. In fact, the satellite eclipses provide 
a simple and effective means of scanning the Jovian atmosphere to determine its physical 
properties as functions of height, latitude, longitude, and time. 


^Researched during the tenure of an NRC-NAS postdoctoral resident research associateship at the Institute for Space 
Studies. 
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CHAPTER 2 


JUPITER AND THE GALILEAN SATELLITE SYSTEM 


A. Physical Parameters: Mean Satellite 

The basic physical parameters of Jupiter and the Galilean satellite system are listed 
in Tables 2. 1 and 2. 2, respectively. These data have been obtained from Allen (1963) 
and Melbourne et al. (1968). 


Table 2. 1 — Physical parameters of Jupiter and its orbit. 


Parameter 

Value 

Equatorial radius of planet 

71,350 km 

Elliptic ity of planet 

0.062 

Inclination of plane of equator to plane of orbit 

3°05’ 

Surface gravity (attractive) 

2600 cm sec -2 

Surface gravity (equator; centrifugal) 

-225 cm sec" 2 

Semimajor axis of orbit 

5. 2028 AU 

Sidereal period 

11. 86223 tropical 

Synodic period 

years 

398. 88 days 

Eccentricity of orbit 

0.048 

Mean daily motion 

299". 1278 

Mean orbital velocity 

13. 06 km/sec 

Inclination of orbit to ecliptic 

1°18' 

Angular diameter of Sun observed from 
orbit of Jupiter (distance from Sun:' 5. 2028 AU) 

368". 9 


To facilitate demonstration of the eclipse technique for inferring the physical prop- 
erties of the Jovian atmosphere, it was decided to use a fictitious mean satellite rather 
than treat each satellite individually. The physical parameters of this satellite (Table 
2. 3) are the averaged parameters of the four Galilean satellites. Note that the mean 
satellite is assumed to revolve around Jupiter in a circular orbit exactly in its equatorial 
plane, to be spherical, and to have uniform reflectivity over its apparent disk. In addi- 
tion, because the difference between the sidereal and synodic periods is negligible, it is 
neglected in the subsequent discussion. For convenience, certain simplifying assumptions 


3 




are made concerning Jupiter itself. The planet is assumed to be spherical and to re- 
volve around the Sun in a circular orbit in the plane of the ecliptic. 


Table 2. 2 — Physical parameters of the Galilean satellites. 



Satellite number 


I 

(Io) 

n 

(Europa) 

m 

(Ganymede) 

IV 

(Callisto) 

Distance from center of 
Jupiter (km) 

4.22 x 10 5 

6.71 x 10 5 

1.070 X 10 6 

1.883 x 10 6 

Sidereal period (days) 

1.769138 

3.551181 

7.154553 

16.689018 

Inclination of plane of orbit to 
Jovian equator (deg) 

0 

0 

0 

0 

Eccentricity of orbit 

0 

.0003 

.0015 

.0075 

Mean orbital velocity (km/sec) 

17.35 

13.71 

10.88 

8.22 

Radius of satellite (km) 

1670 

1460 

2550 

2360 

Apparent visual magnitude, 
m v , at mean opposition 

+5.5 

+5.8 

+5.1 

+6.3 


Table 2.3 — Physical parameters of fictitious satellite. 


Parameter 

Value 

Distance from center of Jupiter (km) 

1.012 x 10 6 

Sidereal period (days) 

6.577 

Inclination of plane of orbit to 
Jovian equator (deg) 

0 

Eccentricity of orbit 

0 

Mean orbital velocity (km/sec) 

11.19 

Radius of satellite (km) 

2010 

Apparent visual magnitude, m v , 
at mean opposition 

+5.54 


B. Eclipses of the Satellites: Basic Geometry 

The geometry of an eclipse of the mean satellite is shown in Figures 2. 1, 2. 2, and 
2.3. Because the ingress light curve is far more easily observed than the egress light 
curve, the discussion will be restricted to the ingress region of the orbit of the satellite. 
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Satellite at Point of Ingress 
to Jovian Shadow 



Figure 2.1— Geometry of an eclipse of the mean satellite. The section shown is that through both the center of 
the Sun and the center of Jupiter, perpendicular to the plane of the Jovian orbit. For simplicity, the special case 
in which the Jovian polar axis NS lies in the plane of the section is shown. In this case, the angle 5 is equal to the 
angle 77, the inclination of the Jovian equator to the plane of the Jovian orbit. In the general case, the angle 8 is 
equal to 77 cos (/— / 0 ), where / is the mean anomaly of the planet, and / Q is the corresponding mean anomaly 
when the north pole of Jupiter is tilted at its greatest angle toward the Sun. The diagram is not to scale. 



Figure 2.2— Geometry of an eclipse of the mean satellite. 
The eclipse shown in Figure 2.1 is viewed perpendicular to 
the plane of the orbit of the satellite. The diagram is not 
to scale. 


To study the effect of refraction by the 
Jovian atmosphere on the ingress light 
curve, it is necessary to be able to calcu- 
late, for any chosen instant during ingress, 
the perpendicular distance p from the cen- 
ter of the satellite to the edge of the Jovian 
shadow. The sharp edge, used for refer- 
ence, would be produced if the planet had 
no atmosphere and the Sun were reduced 
to a point source. With reference to Fig- 
ure 2.3, the coordinates (x, y) of the 
center of the satellite must be defined as 
a function of time t. Time is measured 
from the instant (t = 0) the center of the 
satellite is on the reference edge of the 
shadow at coordinates Xq, y 0 . If r< r s , 
time is counted as positive; if r> r 3 , time 
is counted as negative. The quantity r is 
the distance from the center of the satellite 
to the center of the Jovian shadow when the 
satellite is at a distance from the limb 
of Jupiter. The quantity r s is the radius 
of the point -source geometrical shadow of 
the planet at the same distance from the 
limb. The quantity d 0 is the distance from 
the center of the satellite to the limb of 
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Figure 2.3— Geometry of an eclipse of the mean satellite. The eclipse shown in Figure 2.1 is viewed here from the 
center of Jupiter along the prolonged radius vector from the Sun. The section is taken perpendicular to the axis 
of the shadow at a distance d 0 from the center of Jupiter. The diagram is not to scale. 

Jupiter at ingress, measured along the length of the shadow. To a good approximation 
we can assume that d 0 is equal to r or bit, the radius of the orbit of the satellite. This 
assumption is weakest when T" or bit is small. However, even for Jupiter I the error in 
<3 0 is never greater than 2 percent, and the approximation is, therefore, adequate for our 
discussion. 

The quantity p at any chosen instant is given by 


p = 


r s 


-V * 2 +r 2 , 


(2.1) 


whex’e x , y are the coordinates of the center of the satellite at time t. 
is given by 


1 s 1 Jupiter 


1 + 


'0 


R 


Jupiter 


The quantity r s 


(2.2) 


where /"Jupiter is the radius of Jupiter and ^Jupiter is the distance of Jupiter from the 
Sun. We assume that d 0 is constant for the duration of the ingress light curve. The cor- 
responding uncertainty in the determination of r s never exceeds 5 parts in 10 4 . 

The equation for the reference edge of the Jovian shadow is 


x 2 +y 2 


(2.3) 


and the equation for the orbit of the satellite, projected on the (x , y) plane, is 
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x 2 + 


sin 2 6 orbit ' 


The angle 8 is defined as 


5 = 7? cos ( / ~ / Q ) , 


(2. 4) 


(2.5) 


where 


71 is the inclination of the Jovian equator to the plane of the Jovian orbit 
1 is the mean anomaly of the planet 

_Z 0 is the corresponding mean anomaly when the north pole of Jupiter is tilted at 
its greatest angle toward the Sun. 

The coordinates (x 0 , y 0 ) may be obtained by solving Equations 2. 3 and 2. 4 simultaneously. 

The coordinates ( x , y) at any time t during ingress may be calculated in the follow- 
ing manner. With reference to Figure 2. 2, the angle o; 0 is defined such that 


Q!q = sin 


orbit 


We define another angle, oq , such that 


= “ vt. 


( 2 . 6 ) 


(2.7) 


where v is the mean orbital velocity of the satellite. If we write 

a = + cq , 


it follows that 


X 1 orbit 


sin ct. 


( 2 . 8 ) 


(2.9) 


It then follows from Equation 2. 4 that 


; ; = 



sin 6 . 


( 2 . 10 ) 


Eclipses of the Galilean satellites occur very frequently. For Jupiter I, II, and III, 
an eclipse always occurs during each revolution of the satellite around the planet. In the 
case of Jupiter IV, however, an eclipse does not occur under all circumstances. During 
the course of the Jovian year, the planetocentrie declination of the Sun varies between 
±3°05 ! , depending on how much the polar axis is tilted toward or away from the Sun. The 
planetocentrie declination of the Sun depends on the location of the planet in its orbit 
since the direction of the polar axis is very nearly fixed in space. The large distance of 
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Jupiter IV from the center of the planet prevents the occurrence of eclipses during two 
periods when the planetocentric declination of the Sun is either more northerly or more 
southerly than 2° 10'. In the former case, Jupiter IV passes north of the shadow; in the 
latter case, it passes to the south. Eclipses of the other three satellites continue to oc- 
cur during both these periods because their distances from the center of Jupiter are 
much smaller. 

At all points in the Jovian orbit, eclipses occur so frequently that scarcely a day 
goes by on Earth when one or more does not occur. The eclipses are fairly uniformly 
distributed in time, largely because of a particular relation that exists between the mo- 
tions of Jupiter I, II, and III, This relation is exact and permanent and can be shown 
mathematically to be a necessary consequence of the mutual gravitational perturbations 
between the three satellites. If 0 1# 0 2 , and 0 3 are the mean longitudes of the satellites 
in their orbits, measured from some initial radius through the center of Jupiter, then 

0^ 203 - 302 = 180°. (2.11) 

This relation ensures that all three satellites can never exhibit the same phenomena at 
the same time. Because observations of the eclipses can be made throughout the Jovian 
orbital period, it is possible to monitor the physical structure of its atmosphere quite 
thoroughly in time. Diurnal variation of the atmosphere may also be monitored in a 
crude fashion. The ingress light curves contain information concerning the physical 
structure of the atmosphere at dusk; the egress light curves contain the corresponding 
information at dawn. 

Variation in the physical structure of the atmosphere in both latitude and longitude 
may also be monitored. The atmospheric structure may be scanned quite uniformly in 
longitude since the rotation period of Jupiter is uncorrelated with the orbital periods of 
the satellites. Consequently, each ingress light curve contains information concerning 
the atmosphere at a different longitude. The variation in the planetocentric declination 
of the Sun during the orbital period, caused by the Jovian polar axis not being perpendic- 
ular to the plane of the Jovian orbit, permits the atmospheric structure to be scanned 
completely in latitude twice during the same period. Reference to Figure 2. 1 and Tables 
2. 1 and 2. 2, together with the application of simple geometry, shows that the ranges of 
latitude scanned by the four satellites during the Jovian orbital period are +18° 30' by 
Jupiter I, ±30°18' by Jupiter II, ±53° 37' by Jupiter III, and ±90°00' by Jupiter IV. 

The duration of ingress increases as the absolute value of the planetocentric declina- 
tion of the Sun increases, while the total duration of the eclipse from ingress to egress 
decreases. Reference to Figure 2.3 shows that the satellite strikes the edge of the shad- 
ow at a progressively shallower angle as jc 0 decreases and increases. In particular, 
Jupiter IV can strike the shadow at any angle from normal to tangential. 

To simplify the treatment of the geometry of the eclipse as far as possible in the 
subsequent discussion, it will be assumed that the angle 5 is equal to zero so that at 
ingress the satellite intersects the shadow along the direction of the normal. It follows 
that (Figure 2. 3) the projected satellite orbit will be reduced to a straight line since the 
plane of the orbit then contains the line of sight from Earth. 



CHAPTER 3 


THEORY OF A SATELLITE ECLIPSE: INGRESS LIGHT CURVE 


A. Introductory Remarks 

Before discussing the theory of the ingress light curve in detail, it will be helpful 
to describe briefly the view of the Sun and Jupiter as seen by an observer situated on one 
of the Galilean satellites as it enters the Jovian shadow. If Jupiter had no atmosphere 
above its cloud level, light from the limb of the Sun that was apparently closest to the 
planetary limb would be progressively occulted by the clouds as the satellite gradually 
entered the shadow. At any point during ingress, the Sun would appear as a truncated 
disk. Eventually, the entire disk would be occulted, and no sunlight would reach the 
satellite until it began to leave the shadow at egress. In reality, of course, Jupiter has 
a substantial gaseous atmosphere above its cloud tops that will strongly influence the 
nature of the solar eclipse phenomena seen by an observer on the satellite. First, the 
apparent occultation of the solar disk will be delayed by refraction of the sunlight in the 
atmosphere above the clouds. Second, the shape of the refracted solar disk, and its 
apparent brightness, will be strongly affected by differential refraction. The refracted 
solar disk will appear compressed in a direction perpendicular to the limb of the planet, 
and its brightness will be reduced. Third, both the apparent brightness and color of the 
refracted solar disk will be affected by scattering. If the atmosphere above the clouds 
is clear and gaseous, Rayleigh scattering will both redden the solar disk and reduce its 
intensity. If the atmosphere is not perfectly clear, but contains in addition a haze of 
cloud particles, further extinction of light and modification of its color will result from 
scattering by aerosols. Fourth, throughout the eclipse, the atmosphere of Jupiter above 
the cloud tops will be visible to the observer as a narrow bright ring because of multiple 
scattering of sunlight in the illuminated hemisphere. To the observer, the disappearance 
of the Sun behind the limb of the planet will be quite similar to a sunset on Earth as seen 
from a manned spacecraft in near -Earth orbit. The appearance of the solar disk will be 
the same whether the light rays are prevented from reaching the observer by their in- 
ability to penetrate the thick cloud layer or by critical refraction occurring at some 
height above the cloud tops. 

We will discuss here only the illumination of the satellite by either the unrefracted 
or refracted solar disk. The possible effects on the ingress light curve of haziness in the 
atmosphere above an ill-defined cloud level and the bright atmospheric ring surrounding 
the planet will be discussed later in the paper. For simplicity, we will assume that the 
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atmosphere above the cloud tops is perfectly clear and gaseous, that the cloud tops are 
at a very sharply defined level in the Jovian atmosphere, and that the clouds themselves 
are completely opaque. 

B. Standard Model of Jovian Atmosphere 

A knowledge of the physical properties of the gaseous atmosphere above the cloud 
tops is essential for determining quantitatively the effect of the atmospheric refraction 
on the ingress light curve for the mean satellite. Such information has been obtained 
primarily by spectroscopy. 

Recently, Fink and Belton (1969) analyzed three of the H 2 quadrupole lines in the 
spectrum of Jupiter, viz. the 3-0 S(0), the 3-0 S(l), and the 4-0 S(l) lines. A re- 
flecting layer model was used in their interpretation to obtain both an abundance of H 2 
and a mean temperature for the atmosphere above an effective level of reflection. Their 
analysis gave an H 2 abundance of 67 ±17 km atm above the clouds, and a mean tempera- 
ture of 145° ±20°K. Their work was an extension of that by Owen and Mason (1968) who 
also analyzed two of the H 2 quadrupole lines, viz. the 3-0 S(l) and the 4-0 S(l) lines. 

They derived the abundance of H 2 above the clouds for each line in turn with the aid of 
a reflecting layer model. The results for both lines were consistent. The 3-0 S(l) line 
gave an H 2 abundance of 85 +15 km atm; the 4-0 S(l) line gave 90 ±17 km atm. 

Owen and Mason discussed the problem of the relative abundance of the elements in 
the Jovian atmosphere. Evidence was presented to indicate that the abundances are 
similar to those in the Sun. Molecular hydrogen appears to be the major constituent of 
the atmosphere. An analysis of the collisional broadening of lines in the methane bands 
in the Jovian spectrum shows that the line widths are consistent with an essentially pure 
H 2 atmosphere. There remains, however, a small discrepancy between the theoretical 
and observed line widths, the theoretical width being smaller. This discrepancy may be 
removed by including some helium with H 2 in the collisional broadening theory. The 
maximum amount of helium above the clouds, consistent with the observed line widths, 
is 13 km atm. In that case the H 2 /He ratio is ~7. The corresponding H/He ratio for the 
Sun, given by Biswas and Fichtel (1965), is ~ 11, with a large uncertainty. Methane and 
ammonia, the only gases besides H 2 detected in the Jovian atmosphere, are very minor 
constituents by numbers. 

In application of the theory of the ingress light curve, a standard model will be used 
to describe the atmosphere above the clouds. The effects of variation in the standard 
model will also be studied. The basic assumptions made for each model will be that the 
atmosphere is spherically symmetric, in hydrostatic equilibrium, homogeneous, iso- 
thermal, and composed only of hydrogen and helium. The parameters of the standard 
model are based on the results obtained by Fink and Belton (1969) . The temperature T 
of the atmosphere will be assumed to be 145°K, and the H 2 /He ratio by numbers will be 
taken as 5. 1. The total pressure P 0 at the cloud tops may be derived if the scale height 
of the atmosphere is known. For the scale height, we will adopt the value of 8. 3 km ob- 
tained by Baum and Code (1953) from a stellar occultation. The corresponding value of 
P 0 will be 5. 14 x 10 6 dynes cm -2 . The total particle number density at the surface, 
may be readily calculated from P 0 and T. From the ideal gas law we have 
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p 

TVg = = 2.57 X 10 20 cirf 3 , 

where h is Boltzmann's constant, 

C. Optical Parameters of a Gaseous Atmosphere 


(3. 1) 


Refractivity 

The standard formula for determining the refractivity of a gas at STP is 

( n ~ l)gTP = d ^1 + , (3, 2) 

where X is the wavelength and A and B are constants for the particular gas. Values of A 
and B were tabulated for H 2 and He by Allen (1963). 

The refractivity [n(r)-l] at distance r from the center of the planet is proportional 
to the number density N (r) of the atoms or molecules responsible for the refraction. It 
follows, therefore, that 


[»M- 1] I) 


STP’ 


(3.3) 


where L is Loschmidt's number. 

Specifically, the refractivity at the cloud tops is 


l>0'o) - 0 


^o) 

L 


(n l) sxp . 


(3.4) 


Since we know that 

N(r) = N(r Q ) exp , (3.5) 

where h, the height above the cloud tops, is given by 

h = r-r Q , (3.6) 

r 0 being the distance from the center of the planet to the cloud tops, and H is the scale 
height of the atmosphere, it follows that 

[»(>•)- l] = 0(r Q )- 1] exp . (3. 7) 

Rayleigh Scattering Cross Section 

From Allen (1963), the Rayleigh scattering cross section a may be written as 
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(3.8) 


a =- 


12 87T 5 a 2 


3 A 4 


where a is the polarizability of the atom or molecule. Allen gives a further formula 
which relates the polarizability to the refractivity at STP. We have that 

(n- l) STp = 1.689 X 10 2 V (3.9) 

It is easy to show with the aid of Equation 3. 2 that 



(3.10) 


where a 0 , a*, and 0 % are constants, readily calculable for any atom or molecule from 
the data given by Allen. 


Effective Refractivity and Effective Rayleigh Scattering Cross Section 

If the gaseous atmosphere above the clouds is composed of several different gases 
(e. g. , H 2 and He), it is necessary to use both an effective refractivity and an effective 
Rayleigh scattering cross section. If the total number of atoms and molecules per cubic 
centimeter is written as -^(^)total> then 


^)total = ^)H 2 + ^)He> 


(3.11) 


where N (r) H and /'/(r)j[ e are, respectively, the number densities of H 2 and He at the par- 
ticular point in the atmosphere under consideration. If the effective refractivity is 
written as C^( r )-l]effective> we have 


[«('•) “ 1 ] effective = - 0 H 2 +/ 2 ( r ) ["W “ 1 ] He- (3 ‘ 12) 


The quantities [n(r)-l]H 2 and [n(r)-i] He are, respectively, the refractivities of H 2 and He. 
The refractivities may be calculated with the aid of Equations 3. 2 through 3.7. In each 
equation the number density of atoms or molecules is taken as ■^(rJtotal* The parameters 
/j(r) and / 2 (r) are defined by 


mn 2 

fl(r)= N( r ) H2+ N(r) m 


(3.13) 


and 


/ 2 00 = 


^•) He 

AO') H 2 + A Tr) j j e 


(3.14) 


If the ratio N(r)}i 2 /N(r)iie is assumed to be invariant with height in the atmosphere, / j 
and / 2 will be constant with height also. If the effective Rayleigh scattering cross sec- 
tion is written as o e ff ec tive> we have 
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(3.15) 


^effective A 2 ^ ^2 a He’ 

where CTfj and are, respectively, the Rayleigh scattering cross sections for H 2 
and He. 2 


D. Refraction Through a Planetary Atmosphere 


The effect of refraction on the transmission of light through a planetary atmosphere 
has been discussed by Pannekoek (1903) and Fabry (1929). Both considered only the 
special case of the occultation of a star by a planet as observed from Earth. Attention 
was given to the determination of the angle of refraction, the path length through the at- 
mosphere, and the attenuation by differential refraction. Here the theory of refraction 
in a planetary atmosphere is given in greater detail. The theory has been generalized 
to treat the illumination of a finite satellite by light refracted through a planetary at- 
mosphere from a finite solar source. 


The problem of determining the path of a ray of light through a planetary atmosphere 
reduces to finding the curve T for which the travel time T is a minimum. The travel 
time T is defined by 



ds . . 

— = minimum, 
v 


(3.16) 


where ds is an element of path length, and v is the velocity of light in the atmosphere at 
the corresponding point along the path. The velocity v is given by 


(3.17) 


where n is the refractive index, and a is the velocity of light in a vacuum. The model 
atmosphere adopted for the refraction problem, like the standard model of the Jovian at- 
mosphere, is assumed to be spherically symmetric, homogeneous, in hydrostatic equi- 
librium, and to have a scale height H constant with altitude in the atmosphere. The 
refractive index may, therefore, be written as n(r) since it is a function only of the dis- 
tance from the center of the planet. The element of path length ds may be described by 
polar coordinates as defined in Figure 3.1. We have 


ds 2 = (rdd ) 2 + (dr) 2 . 


(3. 18) 


Equation 3. 16 may, therefore, be rewritten as 


T = 



(3.19) 


where a and b represent the limit in r of the curve F. 
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Figure 3.1— Schematic representation of refraction in the Jovian atmosphere. 


The integral in Equation 3. 19 is of the general form 



b 

®{x,y>y')dx 


(3.20) 


where, in the usual notation, 

- dy 

y = Y x - (3. 2i) 

From the calculus of variations, the Euler -Lagrange differential equation for the ex- 
tremum of J is 


The solution, y = f{x), of Equation 3. 22 gives the equation of the path F 0 in the x, y - 
plane which makes J an extremum. 

Equation 3. 19 is of the particular form 



b 

<b(x,y')dx, 


(3. 23) 


which results when § does not contain y explicitly. In this case the corresponding Euler- 
Lagrange differential equation is simply 


d_ 

dx 



(3.24) 
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It follows that the equation whose solution will determine the path of the refracted light 
ray through the atmosphere is 


4(m)=o, 

drydd J 


(3. 25) 


where 


and 


dr 


f=f(r,d) 


f [Kf) 2+ i 


1/2 


(3. 26) 


(3.27) 


Equation 3. 25 then becomes, after multiplying through by o, 


d j , , T 2 (ddV ^ ,T 

{dT/ +1 J 


1/2 


2 dQ_ I 
dri 


0. 


(3. 28) 


If we integrate Equation 3. 28 once, use to represent the constant of integration, and 
then rearrange the resultant equation, we obtain 


dr r 2/ \ 2 /->2n 1/2 
r \tS(rY - C\\ 

The constant Q x may be evaluated immediately by noting the end conditions correspond- 
ing to a and b , which are the limits in r of the curve. At a, r = Ft, n(r) = n{B), and 
9 = 0. The quantity R is the distance of closest approach of the refracted light ray to 
the center of the planet. If the ray of light is to traverse the atmosphere, we must have 
the condition that 


R>r 0 , (3.30) 

where r 0 is the distance either from the cloud tops or from the level of critical refrac- 
tion, to the center of the planet. If the level of critical refraction is above the clouds, 
the value of r 0 corresponds to that level; vice versa, the value of r 0 corresponds to the 
level of the cloud tops. At b, r = ®, Mr) = 1, 0 = 9^. It follows that 


Cj = n(R)R . 


(3.31) 


The equation to be integrated, therefore, is 

*b 


■i 


n(R)R dr 


' a r[n 2 {r)r 2 - n 2 {R)R 2 } 1 ^ 2 


(3.32) 
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To solve Equation 3.32, a change in variables is made. Let 


it = rn(r ) = rn 
= Rn ( R ) 

du = [ rn'(r) + «(/•)] d/- = (rn + ») o'/', 


where 


n'(r) 


dn(r ) _ dn 
dr dr ' 


Equation 3. 32 may be rewritten as 

0 = 


f 


nu Q du 


\ u(u 2 - + n) 

where u 0 and iq correspond to a and b, respectively. 


(3.33) 


(3 . 34) 


Equation 3. 34 becomes 


f 1 ( | _ m '! n \ 

w 0 / \ 1 + rn'/nj 


u(u 2 - Uq ) 1 ^ 2 


du 


(3.35) 


or 


n du r u x 

~ U °l u(u 2 -u 2 ) 1 l 2 ~ U °J l , 


m'du 


0 u(u 2 - u 2 ) 112 (n+rn') 


(3. 36) 


If iq is replaced by the first term of Equation 3. 36 may be readily integrated. If the 
variable is changed by making the substitution 


2 2 2 
P = U l - U* , 


(3.37) 


the integral takes a standard form. Equation 3. 36 reduces, thei’efore, to 


e = t - 


n _ I rju 

2 ° Ju Q u(u 2 -u 2 ) 1 


du 

1/2 (n + rn) 


(3.38) 


The angle of bending 71 from infinity to the point of closest approach to the center of the 
planet is given by 


'o f 


rridu 


0 ll (“ 2 - “ o ) 1/2 (* + m ') 


(3.39) 
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The total angle of bending e suffered by the light ray in traversing the entire atmosphere 
is, therefore, given by 

£ 

II 

(3.40) 

or 

„ / rn'du 

° Ju 0 u(u 2 - Mq) 1/2 (« + rn') 

(3.41) 

Equation 3. 41 may be rewritten as 


e--2n(W f r 2 "T d \ 

Jr n(r ) \n 2 (r)r 2 - n 2 {R)R 1 ] lU 

(3.42) 

The refractivity, [n(r)~ l], is proportional to the number density N{r) of the refract- 
ing particles. It follows that the refractive index n(r) has the form 

r F ( r-r o) 

n(r)= 1 «0' o ) 1J exp f{ + 1 , 

(3.43) 

since 


T “ r o) 

N(r) = N(r 0 ) exp R 

(3.44) 

Differentiating Equation 3. 43 with respect to r we obtain 


n(r)~ H [»(r 0 ) l] exp // J • 

(3. 45) 

For the particular model planetary atmospheres considered in this paper it is always a 
good approximation to write that 

>i(r) = l 

(3.46) 

and 


n(r) - 1 < < 1 . 

(3.47) 

Substituting Equation 3. 45 into Equation 3. 42 we obtain 


f - Oq ) - l] exp - ° 1 dr 

e = -2 n(R)R / HL ° J L H J 

Jr n(r)\n 2 (r)r 2 - n 2 (R)R 2 y /2 

(3.48) 

Therefore , 


/'A f exp («")* 

H ["°° - 1 eXP wI n{r)^\r) r 2 -n\R)R 2 Y 12 ' 

(3.49) 
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Making use of the approximations in Equations 3. 46 and 3. 47 we obtain 

e.2^-[„( ro )-l]exp( fl )i 

(3. 50) 

To solve the integral in Equation 3. 50 we make the substitution 


t-r- R , 

(3. 51) 

and obtain 


e-2[«(/ 0 ) L H JJ o (f 2 +2Rf) l/2- 

(3. 52) 

If we write 


0 

1 

11 

o 

-s: 

(3. 53) 

where h 0 is the height of the light ray above the cloud tops at the point of closest approach 
to the center of the planet, we obtain 

e= 2 [„( f0 )-l]„exp (-„)! ( , 2 + 2Kf)1/2 - 

(3. 54) 

The integral in Equation 3. 54 is of standard form and may be obtained from Erdelyi 
(1953). Equation 3. 54 reduces, therefore, to 

e = 2 [n(r 0 ) - l]( H ) ex P ( H ) K 0 Q ex P (//) > 

(3. 55) 

where K Q (H /H) is one of the modified Bessel functions. From Watson (1958) we know 
that K 0 (x) has the asymptotic expansion 

expM^W^) -)• 

(3. 56) 

For large values of x we can write 


K 0 (x) exp (x) = 

(3. 57) 

In the case of refraction in the Jovian atmosphere we know that 


f».. 

(3. 58) 
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It follows, therefore, that Equation 3. 55 maybe rewritten as 


e = 2 [n(/ Q ) - 1] | exp ] (3. 59) 

or 

e = (if) ' [ /l(r o) “ 0 ex P (if) ■ (3* 6 °) 

Equation 3. 60 may be rewritten as 

e = e Q exp 

where 

/ 2? r/A 1/2 

e ° = ln 

E. Critical Refraction in a Planetary Atmosphere 

A ray of light being refracted through a planetary atmosphere is moving horizontally, 
at its point of closest approach, to the center of the planet. Critical refraction occurs 
when the refractive bending is sufficient to give the light ray a radius of curvature equal 
to the radius of the surface of constant refractive index that passes through that point in 
the atmosphere. For the case of a spherically symmetric, gaseous, homogeneous, iso- 
thermal planetary atmosphere, it is a simple matter to show from Snell's law of refrac- 
tion that the radius of curvature B 0 of a horizontal ray is given by 




[»('o) ~ ! ] • 


(3. 61) 


(3. 62) 


J 

R 


c 


n(R) - l] M 
n(R) W 


(3.63) 


The level in any particular model atmosphere at which a horizontally moving ray of light 
of a particular wavelength suffers critical refraction can, therefore, be readily deter- 
mined. Rays of light that remain entirely above that level while traversing the atmos- 
phere cannot suffer critical refraction. However, rays of light that penetrate deeper than 
the level, of critical refraction will be trapped by the atmosphere and totally absorbed. 

F. Extinction Along the Path Length of a Ray of Light Traversing a Planetary Atmosphere 

The extinction factor for a ray of light being refracted through a planetary atmos- 
phere is exp {~N^a ) , where a is the effective Rayleigh scattering cross section and Nj 1 
is the number of gaseous particles per cm 2 integrated along the path length through the 
atmosphere . Nt is given by 
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r 


(3.64) 


N t = / N(s) ds. 


From Equation 3. 18 we have 


N r =f N(r) r 2 


+ 1 


1/2 


dr. 


From Equation 3. 44 we have 

n t =2 £ /vo ° )exp [ ,2 (£) +i ] 


1/2 


dr. 


or 


n t = 


“Vo )exp (lf)_( exp +1 ] 


1/2 


dr. 


From Equations 3. 29 and 3. 31 we have 
Nj, = 2N(r 0 ) exp 


which becomes 


n t = 


( r o\ f ( r \( n 2 (R)R 2 

\ H IJ R “ P V H )\„\r)r 2 - o 2 (R)R 2 


» ( . o) ex P (A)j[ ,, ( 


a/2 

+ 11 dr. 


1/2 


dr. 


By rearranging the integral in Equation 3. 69 we obtain 


N t = 2N(r Q ) exp j exp gj 


1/2 


n 2 (R)R 2 


dr. 


1 - 


n 2 (r)r 2 


By making use of the approximations in Equations 3. 46 and 3. 47, 


N t = 2 N(r 0 ) exp ( — 


(d>) r 

v H i [ ' 


Let 


r = R + t . 


(3.65) 


(3. 66) 


(3. 67) 


(3. 68) 


(3. 69) 


(3 . 70) 


(3.71) 


( 3 . 72 ) 
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Substituting Equation 3. 72 in Equation 3. 71 we obtain 


Therefore, 


N t = 2N(r Q ) exp 



(R + t) exp 
(, t 2 + 2Rt) 112 


N r = 2 /V(/' n ) exp 



r°° (R + t ) exp 1 

(it) 

W ) 

Jo (r + 2ft/) 1 /2 


(3. 73) 


(3. 74) 


The integral in Equation 3. 74 may be evaluated readily. It has a standard form and its 
solution may be obtained from Erdelyi (1953). Equation 3. 74 becomes 

N t = 2/V(r 0 ) R ex P K , (|) ex P (§-) , (3. 75) 

where Ky (R /H) is one of the modified Bessel functions. From Watson (1958) we know 
that K\ (x) has the asymptotic expansion 


Since we know that 


exp (x) K ] (x) = 


7 r 

2x) 


1/2 



R_ 

H 


» 1 


for the case of refraction in the Jovian atmosphere, it follows that 


N t = 2N(r () ) R exp 



(3. 76) 


(3. 77) 


(3.78) 


which may be rewritten as 

N t = (2 ttRH) 112 N(r Q ) exp • ( 3 - 79) 

G. Attenuation by Differential Refraction for a Ray of Light Traversing a Planetary Atmosphere 

The effect of differential refraction in the Jovian atmosphere on the ingress light 
curve may be examined with the aid of Figure 3. 2. Differential refraction, resulting 
from the rapid decrease in density with increasing height in the Jovian atmosphere, 
causes an apparent defocusing of light passing through the atmosphere from the solar 
source to an observer situated on the illuminated disk of the satellite. With reference 
to Figure 3.2, the radiation flux per cm 2 from the point source S, at a distance 1 from the 
limb of the planet, has been attenuated by the factor ^because of differential refraction 
in the atmosphere. The factor F is given by 
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Screen Locoted ot 
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Satellite into 
Geometrical Shadow 



Figure 3.2— Schematic representation of the geometry of differential refraction in a 
planetary atmosphere. Point light source assumed. 


where ds is the area of the ring surface that would be produced by the unrefracted light 
beam in the distance 1 if the planet did not possess an atmosphere. The quantity ds' is 
the corresponding area of the ring surface which is produced by the refracted light beam 
in the same distance. By elementary considerations ds is given by 


ds = 2n(L+l)8(L + l ) d8 


and ds 1 is given by 


ds' = 2ir 1 /? / dp , 


so that 


ds ' _ 2tt / 2 /3 dj 3 _ l 2 (3 d$ 

ds 2t r(L + l) 2 8d8 (L + l) 2 8d8 ' 


(3.81) 


(3.82) 


(3. 83) 
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For Equation 3. 83 to be valid, both 6 and p must be small, which is the case of interest 
to us. 

From Figure 3.2, 


5 = ■ (3.84) 

where B is the radius of the planet and 2 is the height above the surface at closest ap- 
proach of the narrow beam of light to the center of the planet. In addition, 

S' -R_LL = e - 5 ; (3.85) 

where e is the angle through which the beam of light is refracted. From Equations 3. 84 
and 3.85 it follows that 


n-BC - (L'-Q 

p l l 


, _ L' S' 
l 


- 5 ' 


or 


Using Equation 3. 84 we obtain 


(3 = 


(R+z) 

Z 


e + 5 . 


p _(/?+;) , {R+z) _ 


By differentiating Equation 3. 88 with respect to 2 , we obtain 

dp _l 1 de _ (L + l) de 

dz 1 L dz IL dz 

By differentiating Equation 3. 84 with respect to 2 , we obtain 

db__ J_ 
dz L 


(3.86) 


(3.87) 


(3.88) 


(3.89) 


(3.90) 


By replacing the derivatives in Equations 3.89 and 3.90 by finite differences, we may sub- 
stitute for dp and db in Equation 3. 83 and obtain 


(R + z) 


IL 


(/ +L)- e 


■][ J 


(L + /) de 
IL dz 


(L + iy 


(R+z) 


(3.91) 
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After some simple algebraic manipulations, we obtain 


F = 



eLI 

C R+z)(L + l ) 



cle IL 
dz (L + I) 


If, for simplicity, we write 


then Equation 3. 92 becomes 


IL 

(/+£)’ 


x l e " 


_ (A+2)_ 

X 1 A_ 


(3. 92) 


(3.93) 


(3. 94) 


By differentiating Equation 3. 60 with respect to height, we can readily show that 


de _ e_ 
dz~ H ’ 


Equation 3. 94 may, therefore, be rewritten as 


F = 



(R + z) 


l + 



(3.95) 


(3.96) 


Equation 3. 96 requires Aq e £ (i? + ^) if A is to be always both positive and finite; i. e. , 
the displacement of the beam on the fictitious screen must be less than the radius of the 
planet. Since in the cases of interest to us Xj e«B, Equation 3. 96 may be used with- 
out difficulty. 

In the case where the point source S is replaced by the Sun, the planet is Jupiter, 
and the screen is replaced by one of the Galilean satellites entering eclipse, Equation 
3. 96 may be rewritten as 


V 

('o + V _ 



(3.97) 


or to a good approximation 


F = 




(3. 98) 


The quantity Xi is then given by 



where d x is the distance of Jupiter from the Sun. To a good approximation 


A 


i 


(3. 100) 


since dy » d Q . 

H. Resolution of the Jovian Atmosphere 

The physical structure of the Jovian atmosphere may be well resolved in both latitude 
and longitude by studying the refraction tails of the ingress light curves of the Galilean 
satellites. For the simplified geometrical situation considered in this paper, each of 
the satellites will enter the Jovian shadow in a direction normal to its edge. Resolution 
of the atmosphere in latitude will be good because the diameter of each of the satellites 
is only a small fraction of the diameter of the planet. The angular resolution of the cir- 
cumference of the Jovian disk will be equal to the ratio of the diameter of the satellite 
to the radius of the planet. In the case of the mean satellite entering eclipse, the reso- 
lution in latitude will, therefore, be equal to 3. 2° . In addition, for the simplified geo- 
metrical situation only, the refraction tail of the light curve will span zero range in 
latitude. Resolution of the atmosphere in longitude will also be good because most of 
the refraction and scattering suffered by a ray of light traversing the atmosphere will 
occur near the point of closest approach to the center of the planet. From Equations 
3. 60 and 3. 79, it is clear that most of the interaction of radiation with the atmosphere 
will occur within a path length of ( 277 BH ) 1 / 2 centered on the point of closest approach of 
the ray to the center of the planet. The resolution of the atmosphere in longitude will, 
therefore, be approximately equal to (2 rtH/B) 1 / 2 radians. Numerically, the resolution 
in longitude will be about 1. 6 degrees, if H is 8.3 km. If the geometrical situation is 
more complex, with the satellite entering eclipse in a direction other than normal to 
the edge of the Jovian shadow, the refraction tail of the light curve will span a finite 
range in latitude, and the resolution in longitude will deteriorate. 

The geometry related to the resolution in height of the Jovian atmosphere at any 
particular point on the refraction tail of one of the ingress light curves is shown in 
Figure 3.3. For simplicity we will consider only the period during ingress when every 
point on the disk of the satellite is illuminated by refracted rays of light from every part 
of the solar disk. The entire refracted solar disk, albeit distorted, is visible to an ob- 
server located at any point on the illuminated hemisphere of the satellite. Ingress has not 
yet progressed to the point where the refracted solar disk would appear to be truncated, 
because the limiting angle of refraction exirnit » for ra y s of light grazing either the 
cloud tops or the level of critical refraction, has not yet been reached for any point on 
the disk of the satellite. 

Consider first the case where the level of critical refraction is below the cloud tops. 
A ray of light traversing the atmosphere will be refracted through the maximum pos- 
sible angle e 0 , when its point of closest approach to the center of the planet is an in- 
finitesimal distance above the clouds. In this case we have 


e limit e 0 ' 


(3. 101) 
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Figure 3.3— Schematic geometry to demonstrate height resolution of the atmosphere at any point on the 
refraction tail of the ingress light curve. The effects of the finite sizes of both the solar source and the satellite on 
the spread in the angle of refraction are shown. 


Consider next the case where the level of critical refraction is above the cloud tops. 
Here, a ray of light traversing the atmosphere will be refracted through the maximum 
possible angle e cr iti ca q, when its point of closest approach to the center of the planet is 
an infinitesimal distance above the level of critical refraction. In this case we have 


€ — 6 
limit critical ' 


(3.102) 


For any chosen model atmosphere, it is a simple matter to determine at what level 
critical refraction occurs for light of any particular wavelength, and to determine 
whether that level is above or below the cloud tops. 


The height h 0 of closest approach to the cloud tops for the ray of light joining the 
centers of the disks of the Sun and the satellite may be obtained from 


e c = e 0 eX P 



(3.103) 


where e 0 is the angle of refraction which is required. For the ray to pass through the 
atmosphere, e 0 must be less than both e 0 and e critical* 
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By differentiating Equation 3. 61 with respect to h Q , we may obtain a relation between 
A#, the angular radius of the satellite viewed from the limb of Jupiter; A0, the angular 
radius of the Sun viewed from the limb of Jupiter; and A h 0 , the spread in height of the 
light rays, which illuminate the disk of the satellite, at their closest approach to the 
center of the planet. We have 


de _ e 
dh Q ~~H’ 


which may be written in finite form as 


Ah 0 

° e c 


(3.104) 


(3.105) 


where 2Ae, the required spread in the angle of refraction, follows from 


Ae = A4> + Ad . 


(3.106) 


The height resolution is, therefore, 2Ah 0 . 

Since A 9 is finite, the beginning of the refraction tail to the ingress light curve oc- 
curs when the center of the disk of the satellite is more than one satellite radius within 
the reference edge of the Jovian shadow. For that particular geometry we have 

e, = A4> , (3.107) 

since reference to Tables 2.1, 2.2, and 2 . 3 shows that 

Ad « A<f> . (3.108) 


It follows, in this case, that 


Ah 0 =H, (3.109) 

and the height resolution, therefore, is 2ff. 

As ingress of the satellite into eclipse progresses further, the required value of 
e 0 rapidly increases and the height resolution improves according to Equation 3. 105. 
The resolution continues to improve as e 0 approaches and h Q gradually de- 

creases. Specifically, a height resolution of one scale height is obtained quite close to 
the beginning of the refraction tail when 


e, 2A4> . 


(3. 110) 


Two special regions of the ingress light curve have not been considered here be- 
cause it is very difficult to make use of them in resolving the structure of the Jovian 
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atmosphere. The first region occurs at the very beginning of ingress when part of the 
disk of the satellite is illuminated by direct sunlight and the other part is illuminated by 
sunlight refracted through the Jovian atmosphere. Here the part of the disk directly 
illuminated dominates in determining the character of the light curve. It is extremely 
difficult, therefore, to extract any information about the refraction properties of the 
Jovian atmosphere from this portion of the light curve. The situation does not improve 
until the area of the disk that is directly illuminated approaches zero. The second re- 
gion occurs at the end of ingress when the disk of the satellite is no longer fully illumi- 
nated by refracted solar radiation and is almost totally eclipsed. Some regions of the 
disk are still receiving radiation from a refracted solar disk, which now appears to be 
increasingly truncated. Other regions no longer receive any refracted solar radiation 
at all because the limiting angle of refraction e ]q m it has been reached for all regions of 
the solar disk. At this point during ingress, variations of the reflectivity over the disk 
of the satellite may arbitrarily and appreciably alter the character of the light curve as 
the illuminated area of the disk becomes smaller and smaller. In addition, accurate 
photometry of this region of the light curve would be extremely difficult since the light 
levels are very low. 

Compared with the effects of the finite sizes of the Sun and the satellite, diffraction 
is of negligible importance in limiting the height resolution of the atmosphere attainable 
by the eclipse technique. If we assume that Jupiter possesses no atmosphere, but has a 
very sharply defined limb, and that the mean satellite is reduced to a point, we can ap- 
ply Fresnel's theory for the case of diffraction by a straight edge to determine the limit 
to height resolution imposed by diffraction. For radiation of wavelength 5500A, the 
diffraction limit is of the order of 15 m. However, since the diffraction limit varies as 
A. 1 ^, diffraction severely limits the effectiveness of any spacecraft flyby radio beam oc- 
cultation experiment to resolve the physical structure of the Jovian atmosphere. Such 
occultation experiments normally utilize radiation transmitted in the S-band (1000 to 
5000 MHz). Consequently, the height resolution attainable using spacecraft flybys of the 
planet will be no better than 5 to 11 km, i. e. , of the same order as the scale height. 

I. Effects of Solar Disk Limb Darkening and Finite Size of Satellite 

As viewed from any one of the Galilean satellites, the angular diameter of the Sun 
is very much smaller than the angular diameter of the Jovian disk. Consequently, the 
curvature of the Jovian atmospheric strata can be neglected in determining the flux of 
radiation incident on the disk of the satellite at any instant during ingress. Also, be- 
cause the linear diameter of the satellite is very much smaller than the linear diameter 
of the Jovian shadow, we can neglect the apparent curvature of its reference edge pro- 
jected onto the disk of the satellite. 

To determine the illumination of any selected point on the disk of the satellite, the 
apparent disk of the Sun is divided into linear strips of equal angular width parallel to 
the limb of Jupiter and perpendicular to the shortest line from the center of the solar 
disk to the reference edge of the Jovian shadow. Each solar strip must be considered 
individually to determine its contribution to the illumination of the chosen point. First, 
it must be decided whether the particular selected strip is visible to an observer situated 
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at the chosen point. To obtain this information, the required value of e must be de- 
termined by simple geometrical considerations. If e < 0, it will not be necessary for light 
from the solar strip to pass through the Jovian atmosphere and be refracted in order to 
reach the observer. The observer, therefore, receives direct and unattenuated radia- 
tion from the solar strip. If e > 0, it will be necessary for the light to pass through the 
atmosphere and suffer refraction and attenuation in the process. In particular, if 
e > the solar strip will be completely invisible to the observer. For the case 

where 0 < e ^ ^limiti the theory described earlier may be used to calculate the total ex- 
tinction suffered by the light on its way to the observer. Summation of the contributions 
from all strips enables the illumination of the observer to be determined. 

In the foregoing discussion it is implied that each solar strip is of equal area and 
brightness. In fact, of course, this is not so, and due consideration must be given to 
the circular shape of the solar disk and also to its limb darkening. Each strip must be 
weighted accordingly. Figure 3. 4 shows the Sun and the geometry associated with the 
limb -darkening problem. 



Figure 3.4-The Sun and the geometry associated 
with the limb-darkening problem. 


To a good approximation the solar 
limb darkening is given by 


1 x (0) = / x (0) (1 - u x + u x cos 6), (3. Ill) 


where I\ (0) is the intensity of radiation of 
wavelength X at position 0 on the solar 
disk (the angle 0 is defined in Figure 3.4), 
The parameter labeled here as u-y is the 
limb -darkening coefficient, tabulated as a 
function of X by Allen (1963). We may 
write 


sin 6 = 



(3.112) 


and 



(3.113) 


where r denotes here the solar radius. By 
transformation, we obtain 


X = 



r 


and 



(3. 114) 


where all units of length are normalized to the solar radius as unity. Equation 3. Ill 
may now be rewritten as 

(3.115) 


/ X (0) = G [i+0'V 1 -(* 2 + y2 )_ 
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where 


G = / x (oxi-« 1 ) 


(3. 116) 


and 

r-Ti \)- < 3 - 117 > 

We can define the quantity A(X) dX to be the fraction of solar radiation in a linear 
strip across the solar disk between X and X + dX. We normalize A(X) dX so that the total 
brightness of the apparent solar disk is unity. In that case we have 

+ 1 

A(X)clX= 1. (3.118) 



We can then write 


f ' f~i + p'\Ji - (x 2 + r 2 )] 

Jy= 0 L J 


A(X) 


l~(X 2 + Y 2 ) I dY 


r f^ 7 y l + ^ l _ (X 2 +Y 2 ) l 

Jx =- 1 Jy = o l j 


(X z + Y l ) \dYdX 


(3.119) 


The evaluation of these integrals is straightforward, and their solution gives 


A(X) 



(3. 120) 


To determine the total brightness of the satellite during ingress, its disk is divided 
into linear strips of equal width perpendicular to the shortest line from the center of the 
satellite to the reference edge of the Jovian shadow. The illumination of a typical point 
within each strip may then be determined as outlined above. If the disk of the satellite 
is uniformly reflective, the contribution of each strip to the total brightness of the satel- 
lite is determined by its area. With reference to Figure 3. 5, the disk of the satellite 
is assumed to be of unit radius and to be divided into n linear strips of equal width. 

The area Bi of a particular strip i is given by 


2yJl-X 2 AX, 


(3.121) 


where the width AX is 2/n, and X i may take any of the values 0, ±2/n, ±4 /n, etc. Be- 
cause the total area of the stips will closely approximate the area of the disk, it follows 
that 

n 

H B i = ' n ( 3 . 122 ) 

i= 1 
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Figure 3.5— The disk of the satellite divided 
into strips of equal width AX. 


if we normalize the area of the disk to 
unity. It follows that 

= AX. (3.123) 

Summation of all contributions permits the 
total brightness of the satellite to be de- 
termined at the particular chosen instant 
during ingress. 

The assumption that the disk of the 
satellite is perfectly circular is an excel- 
lent one. Even when Jupiter is at quadra- 
ture, the defect of illumination is quite 
negligible. The assumption that the disk 
is uniformly reflective is, however, ques- 
tionable. The effect on the ingress light 
curve of an arbitrary distribution of re- 
flectivity over the surface of the satellite 
will be discussed later. 


For convenience in carrying out the 
numerical integrations , both the apparent 
solar disk and the disk of the satellite are divided into 50 linear strips of equal width. 
Numerical experiments, with 200 strips each, showed a negligible improvement in ac- 
curacy in the determination of the ingress light curve. 


Further experiments showed that if less than 20 strips each are used, the accuracy 
of the calculation is decreased. 
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CHAPTER 4 


APPLICATION OF THE THEORY TO THE INFERENCE 
OF THE PHYSICAL PROPERTIES OF 
THE JOVIAN ATMOSPHERE 


A. Ingress Light Curve at Visual Wavelengths: Mean Satellite and Standard Model of Jovian 

Atmosphere 

The ingress light curve for the mean satellite and the standard model of the Jovian 
atmosphere is shown in Figure 4. 1. The wavelength selected for the calculation of the 



Figure 4.1— The ingress light curve for the mean satellite. 
Standard model of Jovian atmosphere: X 5500A. 


light curve was 5500A. In the region of the 
curve -250 sec <t< 250 sec, the bright- 
ness of the satellite is determined largely 
by the contribution from the directly il- 
luminated part of its disk. The refraction 
tail of the light curve does not begin until 
t is approximately 250 sec. For the re- 
gion of the curve with t greater than 250 
sec, the satellite is completely within the 
Jovian shadow. No part of the disk of the 
satellite can then receive solar radiation 
directly. However, for some considerable 
time, the satellite continues to be illum- 
inated by an envelope of light rays from 
the Sun that have been refracted through 
the Jovian atmosphere. As t continues to 
increase, the brightness of the satellite 
steadily decreases. As t advances, each 
of the light rays in the envelope must pene- 
trate increasingly deeper into the atmos- 
phere in order to suffer sufficient refraction 
to reach the disk of the satellite. Deeper 
penetration produces increased attenuation. 
Eventually the limiting angle of refraction 
is reached. Thereafter, as t continues to 
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increase, the light curve rapidly becomes truncated. The location in time of the cutoff 
region in the light curve is directly related to the limiting angle of refraction in the at- 
mosphere. It will be shown later that, for the adopted standard model atmosphere, 
e limit is reached for t > 1500 sec. Consequently, the cutoff in the light curve does not 
appear in Figure 4. 1. Beyond the cutoff region, the satellite is illuminated only by dif- 
fuse solar radiation which has been multiply scattered in the Jovian atmosphere. It 
will be shown later that this source of illumination contributes very little to the bright- 
ness of the satellite. 

The degree of sharpness associated with the edge of the Jovian shadow, for the 
case of the standard atmosphere and a wavelength of 5500A, is indicated by the two curves 
in Figure 4. 2. Curve I refers to the case of a point satellite, with the Sun reduced to a 



Time t (sec) 

Figure 4.2— The degree of sharpness associated with the 
edge of the Jovian shadow. Curve I: Point solar source, 
point satellite. Curve II: Finite solar source, point 
satellite. Standard atmosphere: X 5500A. 


point source. For the purpose of the cal- 
culation, the radius of the satellite was 
set equal to 1 cm, with the angular radius 
of the Sun set equal to 0". 01. Curve II re- 
fers to the case of a point satellite (i. e. , 
radius 1 cm) but with a finite Sun whose 
angular radius is that appropriate for an 
observer at the distance of Jupiter from 
the Sun (i.e. , 184". 45). 

B. Effect of Variation in the Scale Height 
of the Atmosphere 

The effect of variation in the scale 
height of the Jovian atmosphere on the in- 
gress light curve for the mean satellite, 
at a wavelength of 5500A, is shown in Fig- 
ure 4.3. The standard values for both 
the H 2 /He ratio and the pressure at the 
cloud tops were used in the calculation. 
Because the atmosphere is assumed to be 
isothermal and homogeneous, the scale 
height H is independent of height above the 
clouds. The temperature T, associated 
with the chosen scale height, may be ob- 
tained from 

H = — (4.1) 

vg 


where, here, R is the gas constant, g is the surface gravity for Jupiter, and \x is the 
mean molecular weight appropriate for the chosen H 2 /He ratio. The standard value for 
the atmospheric temperature enters into the calculation of each light curve, but only to 
determine the number density N 0 at the cloud tops from the perfect gas law . 
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The light curve that would result if 
Jupiter had no atmosphere is plotted in Fig- 
ure 4. 3 as Curve 1. The other curves, 
numbered 2 through 16, show the refrac- 
tion tails produced by various model at- 
mospheres for a range of scale heights 
between 0 and 100 km. If the standard 
value for the scale height is adopted, the 
corresponding light curve will lie between 
curves 6 and 7. It is clear that the ingress 
light curve is insensitive to the effects of 
atmospheric refraction for extinctions of 
less than 5 or 6 magnitudes. In fact, to be 
able to differentiate clearly between one 
model atmosphere and another, it is nec- 
essary to study the light curve for much 
greater extinctions. It is clear from Fig- 
ure 4. 3 that if accurate observations of the 
light curve could be made for extinctions 
up to and beyond 10 magnitudes , the scale 
height could be very accurately determined, 
provided the H 2 /He ratio by numbers was 
known. Just how such observations can be 
made will be thoroughly discussed later. 

If, however, it is assumed that the refrac- 
tion tail can be located at 10 magnitudes 
extinction with an uncertainty of ±0. 1 mag- 
nitudes, then the accuracy to be expected in the determination of the scale height can be 
estimated. If the true curve lies between curves 6 and 7, the scale height can be deter- 
mined to an accuracy of +1 km if the H 2 /He ratio is known exactly. A precise knowledge 
of the H 2 /He ratio is not strictly necessary, however, because the refraction tail is rel- 
atively insensitive to the helium content of the atmosphere. In fact, as will be shown 
later, if the H 2 /He ratio is near the standard value of 5. 1, an uncertainty in the helium 
content of ±100 percent will produce an uncertainty in the derived scale height of approx- 
imately + 1 km only. In that case, the limiting accuracy attainable in the determination 
of the scale height would be ±2 km. 



Time t{sec) 


Figure 4.3— The effect of variation in the scale height of 
the Jovian atmosphere on the ingress light curve for the 
mean satellite. Standard atmosphere assumed except for 
variation of scale height: X 5500A. 


It may be noted that the scale height computed from the standard values for the tem- 
perature and H 2 /He ratio is inconsistent with the value obtained by Baum and Code (1953). 
If it is assumed that the atmosphere is composed of a mixture of H 2 and He only, with the 
adopted H 2 /He ratio of 5. 1, the mean molecular weight is 2. 33. The corresponding scale 
height for a temperature of 145° K is approximately 21 km, which is substantially differ- 
ent from the value of 8. 3 km (±50 percent) obtained by Baum and Code. Clearly, the 
precise value for the scale height is very uncertain at the present time. Further obser- 
vations of stellar occultations and of eclipses of the Galilean satellites, together with 
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radio wave occultation experiments performed using flyby spacecraft, are required to 
improve our knowledge of the Jovian scale height. 

C. Height Resolution of the Atmosphere 

In Figure 4. 4 the ingress light curve of the mean satellite at a wavelength of 5500A, 
for the case of the standard Jovian atmosphere, is reproduced from Figure 4. 1. Also 

shown is the corresponding curve of height 
resolution of the atmosphere as a function 
of time. The curve of height resolution as 
a function of time begins at the instant dur- 
ing ingress when no part of the disk of the 
satellite is illuminated by direct sunlight. 
As ingress continues beyond this point, the 
height resolution rapidly improves. From 
a value of 25— twice the scale height — 
at the beginning of the refraction tail , it 
improves to approximately 0. 5 H at an 
extinction near 10 magnitudes. Reference 
to the theoretical section on resolution of 
the atmosphere shows that the curve of 
height resolution as a function of time is 
independent of wavelength. 

Figure 4. 5 shows the ingress light 
curve reproduced from Figure 4. 1 together 
with the curve of h 0 as a function of time . 
As before, the quantity fag is the height 
above the clouds of closest approach to the 
center of the planet for a refracted ray of 
light connecting the center of the solar disk 
with the center of the disk of the satellite. 
The curve of fag as a function of time be- 
gins at the same instant as the curve of height resolution as a function of time. Figure 
4. 5 shows that the standard atmosphere can be scanned in height from approximately 
37. 5 km above the clouds downward by studying the refraction tail for radiation of a 
wavelength of 5500A. At an extinction of 10 magnitudes, fag is close to 25 km. The 
value of fag continues to decrease as time increases. The curve of fag as a function of 
time is weakly dependent on wavelength since e is a weak function of wavelength. 

Figure 4. 6 shows the ingress light curve of Figure 4. 1 again, but this time plotted 
with the curve of iY(fag) as a function of time. The quantity N(h 0 ) is the number density 
corresponding to the height fag above the clouds. The curve of N (fag) as a function of 
time is also weakly dependent on wavelength. 

The refraction tail shown in Figure 4. 1 contains information concerning a much 
lower region of the Jovian atmosphere than does the light curve for a stellar occultation. 



Figure 4.4— Height resolution as a function of time during 
ingress of the mean satellite. Standard atmosphere 
assumed. Ingress light curve of Figure 4.1 also plotted. 


36 





Figure 4.5— Height h c of closest approach above the 
clouds of a light ray, connecting the center of the solar 
disk with the center of the disk of the mean satellite, from 
the center of Jupiter as a function of time. Standard 
atmosphere assumed. Ingress light curve of Figure 4.1 also 
plotted. 



Time t(sec) 

Figure. 4.6— Number density N(h c ) corresponding to 
height h c as a function of time. Standard atmosphere 
assumed. Ingress light curve of Figure 4.1 also plotted. 


This refraction tail contains information on the properties of the atmosphere below 
h Q = 37.5 km, corresponding to number densities above 2.83 x 10 18 cm -3 . In the case 
of a stellar occultation observed from Earth, the number density and height, correspond- 
ing to the half -intensity point on the light curve, can be readily calculated by the method 
used by Hunten and McElroy (1968) in their analysis of observations of an occultation of 
Regulus by Venus. The half -intensity point on the light curve observed by Baum and 
Code (1953) for the occultation of a Arietis by Jupiter corresponds to a number density 
in the standard Jovian atmosphere of 1. 32 x 10 13 cm -3 and to a corresponding height of 
139. 3 km. Examination of the light curve obtained by Baum and Code shows that the 
brightness of the star dropped from nearly 100 percent to near zero percent in approx- 
imately 10 sec. Because the apparent velocity of the star normal to the limb of the 
planet at the point of ingress was 12. 0 km/sec, the observed light curve corresponds to 
a range in height of approximately ±60 km centered on a height of 139. 3 km. 


D. Effect of Variation in the H,/He Ratio for the Atmosphere 


Figure 4. 7 shows the effect of variations in the H 2 /He ratio on the ingress light 

o 

curve for the mean satellite at a wavelength of 5500A. Throughout the calculation, the 
other atmospheric parameters were fixed at their standard values. Curves 1 through 6 
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Figure 4.7— The effect of variation in the H 2 /He ratio in 
the Jovian atmosphere on the ingress light curve for the 
mean satellite. Standard atmosphere assumed apart from 
variation in the H 2 /He ratio: A 5500A. 


show the effect of varying the parameter 
#H 2 / H 2 + -^He) from 0 to 1 in steps of 0. 2. 
As before, the quantities N jj 2 and i7jj e re ~ 
fer to the number densities of H 2 and He 
respectively. Because the atmosphere is 
assumed to be homogeneous, the parameter 
#H 2 / (% 2 + l^iie) is invariant with the 
height above the clouds. 

Figure 4. 7 shows that the refraction 
tail of the light curve is insensitive to the 
precise value for the H 2 /He ratio, if the 
helium content of the Jovian atmosphere is 
small. The insensitivity results from the 
relative inefficiency of helium to refract 
and scatter radiation. If the H 2 /He ratio 
is greater than unity it is extremely dif- 
ficult to infer the helium content of the at- 
mosphere, even if the scale height is known 
accurately as a function of height. Con- 
sider, for example, the case where the 
Jovian atmosphere is well represented by 
the standard model. If observations of the 
light curve were made down to an extinc- 
tion of 10 magnitudes, with an accuracy of 
+ 0. 1 magnitudes, the helium abundance de- 
rived from the shape of the refraction tail 
would be uncertain by at least +100 percent 


E. Effect of Variation in the Wavelength 


Figure 4. 8 shows the effect of varying the wavelength on the ingress light curve for 
the mean satellite. The standard model of the Jovian atmosphere was used in the cal- 
culation, Curves 1 through 5 are the light curves for a range in wavelength from 40001 

o 

to 8000A. Each curve has been normalized to eliminate any effect caused by variation of 
the surface reflectivity of the satellite with wavelength. Figure 4. 8 shows clearly that 
the refraction tail is sensitive to variation in the wavelength. The prime reason for the 
sensitivity is the pronounced variation with wavelength of the effective Rayleigh scatter- 
ing cross section. At any selected time on the refraction tail of the light curve, the 
normalized brightness of the satellite is a function of wavelength. The functional de- 
pendence is determined by the optical scattering properties of the chosen model of the 
Jovian atmosphere. 

The degree of sharpness of the cloud tops, as a function of wavelength, is of con- 
siderable interest at the present time. Owen and Mason (1968) have discussed the grow- 
ing body of photometric and spectroscopic evidence that indicates that the cloud tops are 
rather hazily defined. In particular, the height of the effective level of reflection appears 
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to be a function of both wavelength and lati- 
tude. At higher latitudes the cloud level 
appears to be somewhat lower than at lati- 
tudes near the equator. In addition, the 
height of the effective level of reflection 
appears to increase with decreasing wave- 
length. Multicolor photometry of the 
refraction tail will permit the optical scat- 
tering properties of the interface between 
clear and cloudy regions of the Jovian at- 
mosphere to be well studied. The observa- 
tions may be compared readily with 
theoretical predictions for a clear Rayleigh 
scattering atmosphere above a thick cloud 
layer. Disagreement between theory and 
observation is to be expected if additional 
scattering, produced by a haze of cloud 
particles, occurs in this supposedly clear 
region of the atmosphere. The optical 
scattering properties of the cloud particles 
may be inferred and limits placed on their 
probable size and number density as func- 
tions of height. Such information will prove 
useful in checking theories of the nature of 
the clouds in the Jovian atmosphere. 

How high the haze extends above the 
thick cloud level is not known. It is inter- 
esting to speculate that if it should extend to a height corresponding to a number density 
of order 10 13 cm -3 it might severely distort the light curve for a stellar occultation ob- 
served from Earth. During such occultations , differential refraction in the supposedly 
clear regions of the atmosphere high above the cloud layer produces the observed light 
curve. Rayleigh scattering at such levels is minimal. If there should be a significant 
haze of cloud particles in that part of the atmosphere, attenuation by scattering would 
result. Since we would expect the number density of cloud particles in the haze to in- 
crease with decreasing height, the net effect would be to steepen the occultation light 
curve. Unless the effects of such a haze were taken into account in the analysis of the 
light curve, the scale height would be underestimated. Haze may, therefore, account 
for the apparent inconsistency between the value of the scale height obtained by Baum 
and Code (1953) and the value obtained by calculation from the standard values of tem- 
perature and H 2 /He ratio for the Jovian atmosphere. 

F. Effect of Uncertainty in the Radii of the Satellites 

The precise shape of the ingress light curve for each of the Galilean satellites de- 
pends on the value adopted for its radius. How reliably the physical properties of the 



Figure 4.8— The effect of variation in the wavelength on 
the ingress light curve for the mean satellite. Standard 
Jovian atmosphere assumed. Wavelengths between X 
4000A and A 8000A considered. 
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Jovian atmosphere may be inferred from the refraction tails of the light curves is re- 
lated to how accurately the radii of the satellites are known. Although the radii have 
been determined many times by a variety of techniques , the results remain somewhat 
uncertain. 

Camichel (1953) measured the diameters of the satellites by comparing them with 
artificial disks projected in the field of a telescope. The disks were of adjustable size, 
brightness, color, and limb darkening. Dollfus (1954) also measured the diameters with 
a double -image micrometer. The two sets of measurements agree around the following 
values: Jupiter I, 3550 km; Jupiter II, 3100 km; Jupiter III, 5600 km; and Jupiter IV, 

5050 km. 

Interferometric measurements of the diameters of the satellites have been made by 
Michelson (1891a, b), who used a double-slit device in front of the objective of a 12-in. 
refracting telescope. These observations were repeated by Hamy (1899), who used a 
similar technique, but with larger slits to increase the luminosity of the image. Measure- 
ments have also been made by Danjon (1936), who used a half-wave interference mi- 
crometer. Systematic differences appear between the diameters obtained by each 
observer. Diameters measured by Danjon are systematically smaller than those mea- 
sured by Hamy and Michelson; those obtained by Hamy are systematically smaller than 
those obtained by Michelson. Danjon (1944) used an improved reduction technique to 
discuss these interferometric measurements, and succeeded in reconciling them to an 
accuracy of approximately ±5 percent. Basic assumptions in his analysis were that 
each satellite is circular and that its disk is uniformly reflective. His analysis gave 
the following mean diameters for the satellites: Jupiter I, 3300 km; Jupiter II, 2900 km; 
Jupiter in, 4400 km; and Jupiter IV, 3900 km. 

Filar micrometric measurements of the diameters of the satellites have been made 
by many observers. In particular, Barnard (1895) used the Lick 36-in. refractor to 
make an excellent set of measurements. He obtained the following values for the diam- 
eters: Jupiter I, 3946 ±60 km; Jupiter H, 3291 +82 km; Jupiter III, 5726 + 87 km; and 
Jupiter IV, 5383 +56 km. Barnard compared his results with the mean values obtained 
from 9 sets of similar observations made by 9 different observers. His results agree 
closely with those values. Further measurements by Dyson (1895) agree closely with 
the results obtained by Barnard. 

It is by no means clear which of the various techniques for obtaining the diameters 
of the satellites gives the most accurate results. Perhaps the most useful values are 
those obtained by straight averaging of the diameters derived from each of the several 
techniques. The straight mean diameters are: Jupiter I, 3599 +326 km; Jupiter II, 

3097 ±196 km; Jupiter III, 5242 ±733 km; and Jupiter IV, 4778 ±778 km. The error 
quoted is the standard deviation. These errors should give a realistic estimate of the 
uncertainty of the diameters of the satellites. The diameters of satellites I and II ap- 
pear to be known to an accuracy of +5 to 10 percent; the diameter of satellites III and IV 
appear to be known to rather less accuracy, +15 to 20 percent. It will be assumed in 
the subsequent discussion that the diameters are uncertain by up to ±20 percent. 
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Figure 4. 9 shows the ingress light curve for the mean satellite and the standard 
Jovian atmosphere at a wavelength of 5500A for various adopted values for the radius ft 

of the satellite. Curve 4 refers to the case 
where the radius is set equal to 2010 km, 
which is the value listed in Table 2.3 and de- 
noted ftfif for convenience. Curve 1 refers to 
the case of a point satellite, with ft set equal 
to 1 cm. Curves 2 and 6 refer to the cases 
where ft is set equal to 0 . 5 fty[ and 1 . 5 Rtf 
respectively. Curves 3 and 5 refer to the 
cases where ft is set equal to 0.8 7?^ and 1.2 
ftM, respectively. If Rtf should be uncertain 
by up to ±20 percent, the calculated light 
curve would lie between curves 3 and 5. 

The brighter region of the light curve 
is far more sensitive to changes in the 
adopted radius of the satellite than the re- 
fraction tail. This behavior of the light 
curve is a direct consequence of the char- 
acter of the edge of the Jovian shadow, 
depicted by curve 2 in Figure 4. 2, and the 
averaging effect of the finite disk of the 
satellite. Figure 4. 9 shows that the re- 
fraction tail is very insensitive to changes 
in the radius of the satellite for extinctions 
greater than 7 magnitudes. It is clear 
from Figure 4. 9 that the physical proper- 
ties of the Jovian atmosphere can be very 
reliably inferred from the refraction tail, 
even though the radius of the satellite is 
uncertain by up to 20 percent. It is also 
clear that the uncertainty in the radius limits the accuracy required in the photometric 
determination of the shape and location of the refraction tail. If the radius of the satel- 
lite is uncertain by ±20 percent, it is unnecessary to measure the brightness of the satel- 
lite while on the refraction tail to better accuracy than ±0. 1 magnitudes. 

G. Effect of Arbitrary Distributions of Reflectivity Over the Disks of the Satellites 

In reality, we cannot expect the distribution of reflectivity over the disk of a satel- 
lite to be perfectly uniform. We need mention only the case of the Moon, as viewed 
from Earth with the naked eye, to show how much the distribution may depart from 
uniformity. To study the sensitivity of the ingress light curve to the distribution of re- 
flectivity over the disk of a satellite, several artificial and rather extreme situations 
have been studied. Figure 4. 10 shows the ingress light curve for the mean satellite and 
the standard model of the Jovian atmosphere, at a wavelength of 5500A, for various 
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Figure 4.9— The sensitivity of the ingress light curve for 
the mean satellite to the adopted value for its radius. 
Standard Jovian atmosphere assumed: X 5500A. Radius 
varied around the best value R M by up to ±50 percent. 
Light curve for a point satellite also plotted. 
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adopted distributions of reflectivity. 

Curves 1 and 5 both refer to situations 
where one hemisphere of the satellite is 
painted uniformly black (zero reflectivity) , 
while the other is painted uniformly white 
(unit reflectivity). To produce the maxi- 
mum effect on the light curve, we assume 
that the line of demarcation between the 
two hemispheres bisects the disk of the 
satellite and is oriented tangent to the ref- 
erence edge of the Jovian shadow at the 
point of ingress. Curve 1 refers to the 
case where the white section of the disk 
precedes the black section into eclipse. 
Curve 5 refers to the case where the black 
section enters eclipse first. Two other 
less extreme situations have also been con- 
sidered. Their corresponding light curves 
are shown as curves 2 and 4. In both cases 
one hemisphere of the satellite is again 
painted white, but the other hemisphere is 
a uniform grey instead of being perfectly 
black. The reflectivity per unit area of 
the grey hemisphere is one half that for 
the white hemisphere. Curve 2 refers to 
the case where the white section of the disk 
precedes the grey section into eclipse. 

Curve 4 refers to the case where the grey section enters eclipse first. Curve 3 refers 
to the idealized case of a uniformly white distribution of reflectivity over the disk. As 
Harris (1961) first pointed out, the effect on the light curve of changing the distribution 
of reflectivity over the disk of the satellite is similar to the effect produced by changing 
its radius. Consequently, the brighter section of the light curve is again sensitive to 
the physical properties of the satellite. However, the refraction tail is insensitive to 
the distribution of reflectivity over the disk, even for the extremely artificial cases that 
have been considered. Such extreme cases are, of course, very unlikely to occur in na- 
ture. In fact, dark and bright regions are expected to be fairly evenly distributed over 
the surface of the satellite . Such a distribution would cause the average reflectivities of 
both hemispheres to be more nearly equal. Evidence does exist to show that dark and 
bright regions are indeed fairly evenly distributed over the surfaces of the Galilean 
satellites. Maps of their surfaces, based on visual observations made by Lyot, Camichel, 
and Gentili at the Pic du Midi Observatory, have been published by Dollfus (1961). Be- 
sides depicting the coarse surface features of the satellites, the observations showed 
that each satellite's period of rotation was equal to its period of revolution around Jupiter. 
It is fortunate that a detailed knowledge of the distribution of reflectivity over the disks 
of the satellites is not essential. It has been pointed out that the physical properties of 
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Figure 4.10— The effect on the ingress light curve for the 
mean satellite of arbitrary and extreme distributions of 
reflectivity over its disk. Standard Jovian atmosphere 
assumed: A 5500A. 
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the Jovian atmosphere can be reliably studied using photometric observations of the re- 
fraction tails of the ingress light curves accurate to no better than +0. 1 magnitudes. 

Only if the observations are made to greater accuracy need the assumption that the 
disks of the satellites are uniformly reflective be abandoned. 

H. Limits Imposed by Critical Refraction 

If the level at which critical refraction occurs in the Jovian atmosphere is below the 
cloud tops , the location of the cutoff region on the refraction tail of the ingress light 
curve, for each of the Galilean satellites, will be determined by the pressure at the 
cloud tops. If the critical refraction level occurs above the clouds, the refraction tail 
will be prematurely truncated. If the standard model should adequately describe the at- 
mosphere, the height above the clouds, at which horizontal rays of light of wavelength 

o 

5500A are critically refracted, would be 19.2 km, corresponding to a pressure level of 
5. 1 X 10 5 dyne/cm 2 . The eclipse technique for scanning the physical structure of the at- 
mosphere in height could not, therefore, be used to study the region below a height of 
about 20 km above the clouds. There are, however, reasons for believing that the stand- 
ard model may be inappropriate. In particular, the adopted value for the scale height 
may be too small. A larger value would cause a reduction in the adopted pressure at 
the cloud tops. In addition, it would increase the value of the refractivity required to 
produce critical refraction. Both effects would reduce the height above the cloud tops 
at which critical refraction would occur. It is even possible that critical refraction 
might occur within the clouds, in which case the entire region of the atmosphere above 
the clouds could be scanned by the eclipse technique. 

Critical refraction will severely limit the effectiveness of a radio beam occultation 
experiment performed using a spacecraft flyby of the planet. The experiment cannot be 
used to scan the deeper regions of the Jovian atmosphere because the level of critical 
refraction for S-band radiation occurs relatively high in the atmosphere. In fact, critical 
refraction occurs for both S-band and optical radiation at very nearly the same level, be- 
cause refractivity changes only slowly with wavelength. Since refractivity decreases 
with increasing wavelength, the critical refraction level for S-band radiation will be 
slightly lower than the corresponding level for optical radiation. The eclipse technique 
for studying the physical structure of the Jovian atmosphere can, therefore, compete 
very favorably with the radio beam occultation experiment. 
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CHAPTER 5 


OBSERVATIONAL CONSIDERATIONS 


A. Geometrical Problems 

The main geometrical problem is the ability of an observer on Earth to view the 
satellite clearly as it enters the edge of the Jovian shadow. The optimum geometrical 
conditions for observing one of the Galilean satellites at ingress occur when Jupiter is 
at quadrature before opposition. The optimum geometry is shown in Figure 5.1. At 
that time the phase angle of Jupiter will be at its maximum value of approximately 11°, 
and each of the satellites will pass into eclipse at the greatest possible angular distance 
from the Jovian limb. Under the geometrical assumptions made in Chapter 2 of this 
paper the maximum distances from the Jovian limb at which the Galilean satellites enter 
eclipse are, in terms of the apparent radius of Jupiter as unity: Jupiter I, 1. 2; Jupiter 
II, 1.8; Jupiter III, 2. 9; and Jupiter TV, 5. 0. To convert these distances into seconds 
of arc, note that near opposition the apparent angular radius of Jupiter is approximately 
21" of arc. 

Observations of the ingress light curves of the Galilean satellites are possible for 
several months both before and after Jupiter reaches quadrature before opposition. Ob- 
servations may be started, in principle, fairly soon after solar conjunction and several 
months before quadrature. However, certain observational problems occur that are 

related to the geometrical situation. First, the chang- 
ing proximity of the satellite at ingress to the limb of 
the planet plays a role in determining the width in time 
of the useful observational window. Second, the dis- 
tance between Jupiter and the Earth is significantly 
larger before quadrature than after, with the conse- 
quence that the satellite appears fainter before quadra- 
ture and photometric measurements are correspondingly 
more difficult to obtain. Third, the proximity of the 

Figure 5.1— The optimum geometry for the observation from Earth of 
the ingress light curve for the mean satellite. Jupiter is at the quadrature 
point before opposition. Orbits of Earth and Jupiter around the Sun are 
assumed to be coplanar. Plane of orbit of the satellite around Jupiter 
contains the line of sight for an observer on Earth. The diagram is not to 
scale. 
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Sun to the planet also limits the width of the observational window. Observations are 
easiest between quadrature and opposition for several reasons; the angular separation 
of Jupiter and the Sun is larger than between solar conjunction and quadrature; and the 
satellites appear brighter because the distance between Jupiter and the observer has 
steadily decreased. Observations cannot, however, be continued right up to the date of 
opposition because photometric problems become increasingly severe as the satel- 
lites enter eclipse closer to the limb of the planet. To summarize, observations are 
feasible for approximately 3 months per year. The window extends from about 6 weeks 
before quadrature before opposition to about 6 weeks after. 

B. Photometric Problems 

Diffuse Radiation Multiply Scattered from the Jovian Atmosphere 

Mention was made earlier of the illumination of the Galilean satellites, while in the 
Jovian shadow, by the bright ring that appears to surround the planet when viewed from 
within its shadow. The ring is produced by multiple scattering of sunlight in the Jovian 
atmosphere above the cloud tops. During ingress the illumination of the satellite by the 
bright ring is negligible compared with the illumination by the refracted solar disk, as 
will now be shown. 

Consider first the illumination of the mean satellite outside eclipse. Let the solar 
flux at the mean distance of Jupiter from the Sun, and at a visual wavelength of 5500A, 
be Np photons /cm 2 /sec /A. The number of such visual solar photons directly incident 
per second per unit angstrom on the satellite just before ingress will be 

«R 2 s N p , (5.1) 


where B s is the radius of the satellite. 

Consider next the illumination of the mean satellite at the center of the Jovian 
shadow. By the time the satellite reaches this region of the shadow, the refracted solar 
disk will have long since been occulted by the planet and will no longer be visible to an 
observer on the satellite. The satellite will then be illuminated only by the solar photons 
which have been diffusely scattered in the Jovian atmosphere above the clouds. Suppose 
that the atmospheric ring has an effective thickness AZ. The number of solar photons 
incident per second per unit angstrom directly on the ring will be 

2nRAXN, (5.2) 

where Bj is the radius of Jupiter. If it is assumed that the atmospheric ring scatters 
this radiation isotropically, the fraction of these photons which reach the satellite will 
be 


4 nr 2 

S 

or, more simply, 
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( 5 . 4 ) 


R 2 



where r s is the radius of the orbit of the satellite around the planet. 


The effective scattering atmosphere may be assumed to be optically fairly thin, ex- 
tending downward for approximately one scale height from the level at which the optical 
length along a ray path tangent to the limb of the planet reaches unity. Above that level 
the probability of scattering rapidly approaches zero as the number density decreases 
with height. Deeper than one scale height below that level, the probability of a photon 
passing through the atmosphere by multiple scattering is very small because of the rapid 
increase in the optical thickness. For the adopted standard model of the Jovian atmos- 
phere A X is 8.3 km . 


The ratio of the brightness of the satellite outside eclipse to its brightness when 
illuminated only by scattered radiation from the atmospheric ring is 


■nR 2 N 


or, more simply, 


2ttRjAXN p 


2 r‘ 


nR 1 
s 
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(5.5) 


Rj AX ; 


(5.6) 


which is numerically equal to 3. 45 x 10®. It follows that when the satellite is illuminated 
only by the ring it will be 16. 3 magnitudes fainter than outside eclipse. Therefore, il- 
lumination by the atmospheric ring has a negligible effect on the ingress light curves for 
the Galilean satellites for extinctions of less than about 12. 5 magnitudes. 


Accurate Measurement of Satellite Brightness at Low Light Levels 


General Remarks: Introduction to the Area-Scanning Technique 

As mentioned in Chapter 1, the most thorough observational work on the photometry 
of the ingress light curves of the Galilean satellites was carried out many years ago by 
Kuiper and Harris. Their observations have been discussed in detail by Harris (1961). 

It is clear from his discussion that it is extremely difficult to make useful observations 
of the refraction tails of the light curves with conventional photoelectric photometers 
because of scattered light from the nearby bright planetary disk. Kuiper and Harris 
used the 82 -in. -aperture McDonald reflector fitted with a conventional photoelectric 
photometer which had a circular diaphragm of 12" of arc aperture. Because the ap- 
parent angular diameters of the Galilean satellites are approximately 2" of arc, a large 
aperture diaphragm was used to avoid possible vignetting of the image introduced by 
guidance problems, even after the drive rate of the telescope had been corrected for the 
drifts introduced by motion of the planet with respect to the stars, motion of the satellite 
with respect to the planet, and changes in the refraction in the terrestrial atmosphere. 

In principle, a smaller diaphragm, with perhaps a diameter of 4" of arc, could have 
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been used if the periodic errors in the drive screw were less than about 1" of arc. 
Diaphragms with diameters less than about 4" of arc would not have been feasible be- 
cause of variation in the diameter of the seeing disk of up to about 1" of arc. Kuiper 
and Harris took no special precautions to minimize scattered light from the nearby 
planet in the optical system of the telescope. Consequently, they found that even if each 
of the satellites was observed near quadrature before opposition, accurate measure- 
ments of their ingress light curves could not be made fainter than an extinction of about 
6 or 7 magnitudes. In addition, they found that the satellites became undetectable against 
the bright sky background for an extinction of about 9 magnitudes. Below that level the 
signal from the satellite disappeared in the noise of the bright background. The observa- 
tions were made through awide-band yellow filter; Kuiper's earlier visual observations 
had shown that the eye could follow the satellite to about the same limit. The limiting 
magnitude for the visual observations against the bright background was roughly es- 
timated as 14. 5. 


The experience of Kuiper and Harris suggests that even a relatively minor improve- 
ment in photometric technique should allow accurate photometric observations of the 
refraction tails of the light curves to be made. The problem is how to detect and ac- 
curately measure a faint signal against a bright background. Because conventional 
astronomical photometers are not designed for this type of measurement, a new tech- 
nique must be used. A partial breakthrough was achieved by Rakos (1965, 1966) who 
.successfully developed an area -scanning photometer. The advantage of area scanning 
is that, by integrating many scans, the faint signal appears above the noise of the bright 
background. As the integration proceeds, the noise level of the sky decreases while 
the faint signal is accentuated. Very soon the signal appears well above the noise level 
and may be accurately measured. Rakos developed the scanning technique specifically 
to observe eclipses of Phobos by Mars. His objective was to infer the physical proper- 
ties of the Martian atmosphere from the ingress light curve. The observational prob- 
lem he encountered was much more difficult than our own. He had to make observations 
of Phobos when it was only 4” of arc from the bright limb of Mars and when the magni- 
tude difference between planet and satellite, before ingress, was 13. In the case of the 
Galilean satellites entering eclipse, observations maybe made some tens of seconds of 
arc from the bright limb at quadrature before opposition. At such great distances from 
the bright limb, the sky brightness is much less than it is when close to the planet. In 
addition, the magnitude difference between planet and satellite, before ingress, is only 
about 8. The scattered light problem is, therefore, far less severe. 

In the operation of the Rakos photometer, a small region of the sky is optically 
scanned, by means of a rocking mirror technique, through a narrow slit situated in 
front of a Fabry lens/photomultiplier system. In this manner a rectangular strip of 
sky, with dimensions of perhaps 4" x 10", depending on the scale in seconds of arc per 
millimeter at the slit of the photometer, may be repeatedly scanned. Appreciable changes 
in the seeing do not normally occur in a clear sky over such a small area in less than 
~ 0. 1 sec, which is the time required to complete each scan. Individual scans are dis- 
played on an oscilloscope and then recorded on film. Averages are then taken of groups 
of 100 consecutive scans, to reduce the statistical fluctuations of the individual scans. 
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However, because the averaging must be done by hand, it is unavoidably subjective. 

The area-scanning photometric technique has been developed further since the pioneer- 
ing work done by Rakos. Very recently Hall (1968) constructed a twin-beam area- 
scanning polarimeter. He greatly improved the efficiency and sensitivity of the 
area -scanning technique by automatically carrying out the integration of the individual 
scans at the telescope, while the observations were actually in progress. Each of the 
beams of light in the polarimeter was fed to a separate photomultiplier and amplifier. 

The photometric data were then fed into two multichannel analyzers that integrated the 
individual scans. The averaged scan obtained during each integration period was punehed 
on tape and later reduced by computer. The Hall photometer is ideally suited to the 
measurement of the refraction tails of the ingress light curves of the Galilean satellites. 
It will be shown that the necessary observations may be made with this photometer at- 
tached to a telescope of only moderate size, e.g. , only 72-in. aperture. 

Observational Technique 

Ideally, the telescope should accurately track the satellite as it enters the Jovian 
shadow. Allowances should be made for the motion of the satellite with respect to 
Jupiter, in addition to the motion of Jupiter with respect to the stars and the drift pro- 
duced by variation in the refraction in the terrestrial atmosphere during the course of 
the observations. However, it may not be possible to guide the telescope on the satel- 
lite itself when ingress has progressed to the point where the satellite's disk is being 
illuminated only by refracted solar radiation. The satellite may then be too faint to be 
seen easily. The telescope should, therefore, be guided by an offset technique on the 
limb of Jupiter from the very start of ingress. 

The rectangular area of the sky that is scanned by the photometer should be oriented 
in one of several ways, depending on whether the telescope is tracking on the satellite 
itself or on Jupiter. In the ideal case, when the telescope is accurately tracking on the 
satellite, the scanned area should be radially oriented to the planet because the bright- 
ness of the sky background is a function only of the radial distance from the center of 
the planetary disk. Radial orientation would permit the scattered light to be very ac- 
curately subtracted from the signal from the satellite itself. When the telescope is 
guided on the limb of Jupiter, the scanned area should be oriented so that the apparent 
direction of motion of the satellite, with respect to the planet, is exactly along the great- 
est length of the rectangle. 

Care must be exercised in choosing both the dimensions of the scamied area and the 
width of the slit of the scanner. The area scanned should be about 4" of arc wide and 
about 10" of arc long. The width of the rectangle should be kept as small as possible, 
to keep the scattered light level to a minimum. A slit length of 4” of arc will allow for 
fluctuations in the seeing of an order of 1" of arc or less, and will permit small errors 
to be made in the orientation of the scanned area. The slit should be no wider than about 
0". 2, to permit the disk of the satellite to be adequately resolved. Resolution of the 
satellite is desired to maximize the signal/sky ratio. 
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A time resolution of between 10 and 20 sec should be selected so that the precise 
shape of the ingress light curve can. be obtained for each of the satellites. During the 
course of the integration period, it is desirable that the satellite not change its position 
in the scanned area by any appreciable amount. If it does, the reason for making the in- 
tegration is lost. When the telescope is accurately tracking the satellite there is no 
problem. When the telescope is tracking on the limb of the planet, the motion of the 
satellite with respect to the scanned area in a 10 -sec period will typically be less than 
0". 1, a very small fraction of the apparent diameter of the satellite. 

Use of the Hall photometer will enable all these observational problems to be over- 
come very simply. The photometer is fitted with an offset guider. It will be a simple 
matter, therefore, to guide the telescope on the limb of the planet. At ingress, the drift 
direction of the satellite, with respect to the planet, will be closely parallel to the Jovian 
equator. Orientation of the scanned area will, therefore, be simple because the satellite 
will drift parallel to the dark belts visible on the disk of the planet. The scanned area 
can be easily located so that during ingress the satellite moves along the length of the 
rectangle from one end to the other. To orient the scanned area, the whole photometer 
may be easily and accurately rotated around the optical axis of the telescope. 

Photon Statistics and Telescope Aperture 

In the application if an area-scanning photometer to the accurate measurement of the 
refraction tails of the ingress light curves of the Galilean satellites, the length of the 
area scanned is divided into a large number of channels, perhaps 100 or more. As the 
area is scanned, photoelectron counts are recorded in each of these channels and inte- 
grated over many scans, using a multichannel analyzer. A typical record from the 
multichannel analyzer at the end of an integration period is shown in Figure 5. 2. As 
indicated earlier, the signal is accentuated relative to the noise of the sky background. 
The noise of the background is reduced as the number of counts increases. It is im- 
portant to understand that the noise results merely from statistical fluctuations in the 
relatively small number of photons counted. If both the sky transparency and seeing 
conditions are good, it can be assumed that the brightness of a particular point in the 
sky is constant to within 1 percent for a period of up to an hour or so. 

With reference to Figure 5. 2, the area under the signal spike and above the sky 
background is linearly related to the brightness of the satellite. It is necessary to 
know how accurately that area can be determined so that the uncertainty of the measure- 
ment of the satellite brightness may be reliably estimated. The area is proportional to 
the total number of counts recorded above the sky background, when photons from the 
disk of the satellite are passing through the slit of the area scanner. To discuss the 
statistics of photon counting, certain parameters must be defined. Let the length of the 
slit of the area scanner be Aw; the width of the slit of the area seamier be A 1 ; the total 
length of the area of sky scanned be 1 ; the duration in time of each scan be At; the total 
number of scans in each integration period be N; the integration period be T; the flux of 
photons per cm 2 per sec, integrated over the bandwidth of the chosen filter and arriving 
at the top of the atmosphere from, the satellite, be A sa y; the corresponding flux of photons, 


50 



INTEGRATED NUMBER OF COUNTS PER CHANNEL 



Figure 5.2— A typical record from the multichannel 
analyzer at the end of an integration period. 

scan, when light from an equivalent width d 
slit of the area scanner, be n- s ky. 


apparently from each square second of ai*c of 
sky near the satellite, be# s py;the angular 
diameter of the satellite be d sec of arc; 
the ratio of the brightness of the satellite 
to the brightness of a neighboring region 
of sky of equal width d , but of area (d. A w ) , 
be B ; the efficiency of the observational 
system in converting photons incident on 
the top of the atmosphere into photoelec- 
trons in the photomultiplier (counts) be /; 
the diameter of the aperture of the tele- 
scope be D ; the number of photoelectrons 
counted per scan, when light from the 
satellite is passing through the slit of the 
area scanner, be (7Tgj_g + rigpy) ; and the 
number of photoelectrons counted per 
of neighboring sky is passing through the 


For a given set of equipment the accuracy of the measurements is dependent pri- 
marily on the selected period of integration. The relation between T, N, and A t is 
given by 


T = N At . 


(5. 7) 


Clearly, we may write 


«sig =7V sat 



and 

UyVAWrfAH.) /fO//)' 
The basic signal /sky ratio, B, is given by 
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(5.8) 


(5. 9) 


(5. 10) 


Obviously to optimize the signal/sky ratio, the area of sky being sampled simultaneously 
with the satellite must be as nearly equal to the area of its disk as possible. Since Aw 
has been selected to be close to 2 d to ease the guidance problem during the observations, 
the area of sky being sampled is rather more than twice the area of the satellite. 

There are three basic statistical problems involved in the accurate determination of 
the brightness of the satellite. The first is concerned with the accurate measurement 
of the integrated signal from the satellite alone. To study this problem, it can be 
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assumed for simplicity that the brightness of the sky background is zero. After N scans, 
the integrated count for the satellite is By elementary statistics the uncertainty 

in the count is ' \sig. The percentage of accuracy in the determination of the bright- 

ness of the satellite will be 100 /\Z^sig- 

The second problem is concerned with the accurate measurement of the integrated 
signal from the satellite against the sky background. For simplicity, it can be assumed 
that the integrated count for the satellite alone, ilfogig, is known with infinite accuracy. 
After II scans, the integrated count for the sky ^ft s ky, will be known with an accuracy of 

A V^sky- We may define the parameter to be the ratio of the integrated signal from 
the satellite alone to the noise level of the integrated sky count. Therefore, 

Nil ■ 

- -_jg (5.11) 

or 

R i= V A '"sk^ R - (5.12) 

The third problem is to determine the minimum brightness of the satellite when it 
is still just detectable against the noise of the sky background. The parameter 7? 2 can 
be defined as the minimum value of R that is detectable with the particular observational 
system being used. For the signal from the satellite to appear above the noise level of 
the sky after integration over N scans , it is required that 




(5. 13) 


or 


It follows that 


R > 
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Nn 


sky 


(5. 14) 
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(5. 15) 


To demonstrate the feasibility of accurately measuring the ingress light curves to 
low light levels at visual wavelengths with a telescope of only moderate size, the use of 
the Hall twin-beam photometer attached to a 72-in. Cassegrain reflecting telescope 
will be considered. Assume that the telescope is located at Flagstaff, Arizona. Ob- 
servations of the light curves, in two colors simultaneously, may be made provided 
certain minor modifications of the photometer are carried out. In particular, it is 
necessary to restrict the length of the area scanned by changing an appropriate cam in 
the photometer. To facilitate detection of the satellite against the sky background, the 
integration period should be 20 seconds to allow 20 scans to be made at the rate of one 
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per second. A longer integration period will distort the observed shapes of the ingress 
light curves. If the refraction tails of the light curves are to be well defined, accurate 
photometric measurements of each satellite must be made when it has suffered up to 10 
magnitudes extinction. At that level of extinction, it is necessary to determine the 
brightness of the satellite to an accuracy of 10 percent or better. 

It will be assumed that observations are made using the V-filter, which corresponds 
to an effective wavelength of 5500A, and that the dimensions of the area scanned are 
restricted to 4" x 10". The collecting area of the primary mirror of the telescope 
{D = 72 in. ) is 2. 62 X 10 4 cm 2 . The efficiency of the observational system is the product 
of the efficiencies of its individual components, which are as follows: 

(a) Typical transmission efficiency of clear atmosphere =0.7 (altitude of Flagstaff, 
2210 meters; elevation of Jupiter > 30° above horizon, i.e. , atmospheric path < 2 air 
masses). The transmission efficiency is derived from Allen (1963) who gives "V- 
absorption for a clear atmosphere = 0.20 mag/air mass." 

(b) Efficiency of Cassegrain telescope (2 mirrors; aluminum coatings) = (0. 8) 2 
= 0.64. 

(c) Efficiency of Hall photometer (each beam) : 

4 optical components = (0. 8) 4 = 0. 41, 

Wollaston prism beam-splitter = 0. 5, 

V-filter = 0.5, 

Photocathode = 0.1. 

It follows that the efficiency of the observational system is 0.7 X0.64 X0.41 x0.5 X0.1 
X 0. 5 = 0.0046« 0. 5%. 

Consider next the fluxes of satellite and sky photons within the V-filter bandwidth 
which strike the top of the atmosphere before and during ingress. From Table 2.3, the 
apparent magnitude of the mean satellite just before ingress will be My = +5. 54. The 
apparent magnitude at 10 magnitudes extinction will be Wy = +15. 54. From the experi- 
ences of Kuiper, Harris, and the author, in attempting to carry out photoelectric pho- 
tometry of the refraction tails of the ingress light curves when Jupiter nears quadrature 
before opposition, the brightness of an area of background sky equal in area to the disk 
of the satellite, typically should correspond to an apparent magnitude m v = +14. 5. Of 
course, if the aluminum surfaces of the mirrors are badly oxidized, dirty, or scratched, 
the sky brightness will be considerably higher. The flux of radiation, at visual wave- 
lengths outside the Earth's atmosphere, which corresponds to a given apparent stellar 
magnitude m v , may be obtained from the formula given by Allen (1963) 


log 10 / x = *0.4m ); - 8.42 , (5.16) 

where /\is the flux in ergs/cm 2 /sec/A. For m v = +5. 54, log 10 j\ = -10. 62, and = 2.4 
X 10 -11 ergs /cm 2 /sec /A. Let us assume that the equivalent width of the V-filter is 500A. 
The flux F v integrated over the width of the V-filter is, therefore, equal to 1.2 X 10 “ 8 
ergs/cm 2 /sec. The energy of a photon of wavelength 5500A is equal to 3. 6 x 10~ 12 ergs. 


53 



If N v is the number of photons striking the top of the atmosphere in the V -filter band- 
width /cm 2 /sec, then 

/V r = 3.3 X 10 3 photons/cm 2 /sec, for in y = +5.54 , 

TV, = 0.33 photons/cm 2 /sec, for m y = +15.54 , (5„ 17) 

and 

N y = 0.83 photons/cm 2 /sec, for m v = +14.5 . 

For the mean satellite at an extinction of 10 magnitudes, it follows from Equation 5. 8 
that 

n. = 0.33 X 2.62 X 10 4 X 0.005 X (1 X 2/10) (5. 18) 

or 

n • = 8.6 counts/scan, (5. 19) 

for a scan rate of one per second. Similarly, 

"sky = (2 x 0.83) x 2 - 62 X 10 4 X 0.005 X (1 X 2/10) (5. 20) 

or 


« sky = 43- 2 ccmnts/scan. (5.21) 

In Equation 5. 20 the number of photons corresponding to the sky background has been in- 
creased by a factor of 2, because the area scanned has a width twice the diameter of the 
satellite. 

From Equation 5. 19, the accuracy of the measurement of the signal from the satel- 
lite, at 10 magnitudes extinction, will be 13. 1 percent for an integration period of 20 
seconds. From Equation 5. 10 the parameter Rtor 10 magnitudes extinction is 0. 2. It 
follows from Equations 5.12 and 5.21 that R^ i- s equal to 5.9. The accuracy of the meas- 
urement of the integrated signal, from the satellite against the sky background, will be 
about 17 percent. The apparent magnitude of the satellite, when it becomes no longer 
detectable, may be estimated from Equations 5. 15 and 5. 21. Thus, R% is equal to 
0.034, which corresponds to an apparent magnitude for the satellite of about 18, i. e. , 

12. 5 magnitudes extinction. 

It should be borne in mind that the estimates of background sky brightness, efficiency 
of the observational system, and bandwidth of the filter, have been deliberately chosen 
to be conservative. It can be concluded, therefore, that valuable observations of the re- 
fraction tails of the ingress light curves of the Galilean satellites can be made with area- 
scanning photometers of the type developed by Hall attached to telescopes of only mod- 
erate size. 
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CHAPTER 6 


CONCLUDING REMARKS 


A radio beam occultation experiment performed during a spacecraft flyby mission 
to Jupiter is not the most efficient technique available for studying the physical struc- 
ture of the Jovian atmosphere. This paper has described in detail how eclipses of the 
Galilean satellites may be utilized to give a simple, inexpensive, yet extremely power- 
ful method of thoroughly monitoring the physical properties of the Jovian atmosphere 
above the visible cloud layer. The relevant observational data may be obtained with 
telescopes of only moderate size coupled with relatively simple photometers. It appears 
that much of the data required for a complete understanding of the physical structure of 
the Jovian atmosphere may be obtained without the use of spacecraft, but from available 
ground-based astronomical facilities. 


Goddard Space Flight Center- 
National Aeronautics and Space Administration 
Greenbelt, Maryland, December 15, 1969 
039-52-01-01-51 
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